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ABSTRACT 

We use photometric and spectroscopic observations of the echpsing binary 
V69-47 Tuc to derive the masses, radii, and luminosities of the component stars. 
Based on measured systemic velocity, distance, and proper motion, the system 
is a member of the globular cluster 47 Tuc. The system has an orbital pe- 
riod of 29.5 d and the orbit is slightly eccentric with e = 0.056. We obtain 
Mp = 0.8762 ± 0.0048 M0, Rp = 1.3148 ± 0.0051 i?©, Lp = 1.94 ± 0.21 
for the primary and M, = 0.8588 ± 0.0060 M©, = 1.1616 ± 0.0062 i?©, 
Ls = 1.53 ± O.17L0 for the secondary. These components of V69 are the 
first Population II stars with masses and radii derived directly and with an 
accuracy of better than 1%. We measure an apparent distance modulus of 
(m — M)v = 13.35 ± 0.08 to V69. We compare the absolute parameters of V69 
with five sets of stellar evolution models and estimate the age of V69 using mass- 
luminosity- age, mass-radius- age, and turnoff mass - age relations. The masses, 
radii, and luminosities of the component stars are determined well enough that 
the measurement of ages is dominated by systematic differences between the 
evolutionary models, in particular, the adopted helium abundance. By com- 
paring the observations to Dartmouth model isochrones we estimate the age of 
V69 to be 11.25±0.21(random)±0.85(systematic) Gyr assuming [Fe/H] = -0.70, 

^Carnegie Observatories, 813 Santa Barbara St., Pasadena, CA 91101-1292; 
(ian,burley ) @obs . carnegiescience . edu 

•^Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland; (jka,pych)@camk.edu.pl 

"^David Dunlap Observatory, Department of Astronomy and Astrophysics, University of Toronto, P.O. 
Box 360, Richmond Hill, ON L4C 4Y6, Canada; rucinski@astro.utoronto.ca 

^Las Campanas Observatory, Casilla 601, La Serena, Chile; wojtek(S)lco.cl 

^Department of Physics and Astronomy, University of Victoria, P.O. Box 3055, Victoria, B.C., V8W 3P6, 
Canada; dotter@uvic.ca 



-2- 



[a/Fe] = 0.4, and Y = 0.255. The determination of the distance to V69, and 
hence to 47 Tuc, can be further improved when infrared echpse photometry is 
obtained for the variable. 

Subject headings: binaries: close - binaries: spectroscopic - stars: individual 
V69 47 Tuc - globular clusters: individual (47 Tuc) 



INTRODUCTION 



Detached eclipsing double-line binary stars are the fundamental astrophysical laboratory 
for the determination of stellar parameters of mass and radius. Luminosities can be derived 
using measured parallaxes or from empirical color - effective temperature relations. These 
data are the fundamental tests of s tellar evolution models. Many field Population I systems 
are known at solar mass and larger (j Andersenll 1 99 ll ) and modern high accuracy measurements 
of masses, luminosities, and radii o f the component stars are in general agreem ent with 
evolution models (see, for example, iLacy et al.l l2005l . l2008l : IClausen et al.ll2008l ). Similar 
results are obtained for studies of individual binaries in the old open cl usters NGC 188 



flMeibom et al.l 120091 ). NGC 2243 flKaluznv et al.l 120061 ). and NGC 6791 flCrundahl et al. 



20081 ) . Recent wide field photometric surveys have identified numerous low mass systems, 
and for these K and M stars the common theme of a comparison of component properties with 
evolution models is that the models systematically underestimate the ra dii of the components 



i n the s e binaries. Suniniaries of re cent measurements can be found in iLopez-Morales et al. 
j2006h . lL6pez-Moralesl J2007h . and lsiake et~aD J2008h . 



The situation is even less clear for Popul ation II star s. With the exception of CM Dra, 
for which the component masses are ~0.2 Mq (|Lacylll977l ). and the u Cen binary OGLEGC- 
yi7, for which th e analysis is comprornised b y an uncertain determination of the metallicity 
(iThompson et al.ll200ll : iKaluzny et al.l |2002| ) . there are no know n Populatiori II detached 



double-line eclipsing binaries with main-sequence components. iTorres et al.l (|2002| ) used 
interferometric observations of HD 195987 ([Fe/H] ~ -0.6) to derive an orbit and measure the 
masses of the components. The radiative properties of the two components agree with a suite 
of models given some slight modifications of input parameters. While direct measurements 
of the radii of the components are not possible because HD 195987 is not an eclipsing binary. 



^This paper includes data gathered with the 6.5-meter Magellan Baade and Clay Telescopes and the 2.5- 
meter du Pont Telescope located at Las Campanas Observatory, Chile. It is based in part on data obtained 
at the South African Astronomical Observatory. 
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estimates of the radii can be derived from the orbital parallax, the bolometric flux, and the 
estimated effective te mperatures. Her e again the measured radii are larger than the models 
by some ten per cent. iBoyajian et all (120081 ) have measured the angular diameter of the G 
sub dwarf fi Gas ([Fe/H ~ -0.8). This measurement provides a radius when combined with 
the Hipparchos parallax for /i Gas. For this star the models underpredict this measured 
radius by about five per cent. The masses of the components of n Gas are only known to 
about ten per cent, so the model comparisons here are not well constrained. 

There is a clear need to locate and study Population II detached eclipsing binary stars 
to obtain accurate masses and radii of their component stars. Stellar evolution models 
are becoming increasingly sophisticated and are used to fit observed cluster color-magnitude 
diagrams (GMD). The cluster GMD's are themselves improving in quality, wi th homogeneous 



surve ys being conducted with the Hubble Space Telescope (see, for example, ISarajdeni et al. 
20071 ). Gareful empirical tests of these models are an essential next step. 



This is the first paper in a series devoted to the study of detached eclipsing double-line 
binaries (DEB) in Galactic globular clusters with components on the cluster main sequence 
or subgiant branch. The Gluster AgeS Experiment (GASE) has the goal of determining 
the basic stellar parameters (masses, luminosities, and radii) of the components of cluster 
binaries to a precision of better than 1% in order to measure cluster ages and distances, 
and to test stellar evolution models. The methods and assumptions utilize basic and simple 
approache s offere d by the field of eclips ir ig dou ble-line spectroscopic binaries as described in 
Paczyhskil (119971 ) and [Thompson et al.l (120011). Pr evious GASE papers have discussed blue 



straggler systems in u Gen ( Kaluznv et al 
an SBl binary in NGG 6397 lKaluzny et al. 



2007a|) and 47 Tuc (IKaluznv et al.ll2007bf ). and 



(2008) 



The eclipsing binary V69-47 Tuc (hereinafter V69) was discovered by IWeldrake et al. 



(120041 ) during a survey for variable stars in the field of the globular cluster 47 Tuc. They 
presented an /-band light curve for the variable and proposed an orbital period of P = 
5.229 d. The light curve phased with this period shows only one eclipse. The variable is 
located at the top of main sequence in the cluster color- magnitude diagram, and thus is of 
potential great interest for measurements of the cluster age and distance. 

In this paper we report the results of photometric and spectroscopic observations aimed 
at a determination of the absolute parameters of the components of V69. Section [2] describes 
the photometry of the variable and the determination of an orbital ephemeris. Section [3] 
presents the radial velocity observations. The combined photometric and spectroscopic ele- 
ment solutions are given in SectionlHwhile the membership in 47 Tuc is discussed in Sectional 
In Section [6] we compare the properties of the components of V69 to a selection of stellar 
evolution models with an emphasis on estimating the age of the system. Finally in Section [7] 
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we summarize our findings. 



PHOTOMETRIC OBSERVATIONS 



The bulk of tlie pliotometric data were obtained witfi tfie 1.0-m Swope telescope at tfie 
Las Campanas Observatory using tlie 2048 x 3150 pixel SITES CCD camera witli a scale 
of 0.435 arcsec/pixel. These observations were collected during the 2004 - 2007 observing 
seasons. The same set of BV filters was used for all observations. Exposure times ranged 
from 100 s to 360 s for the V filter (average exposure was 120 s) and from 160 s to 360 s 
for the B filter (average exposure was 170 s). An eclipse event was detected on the first 
of a total of 8 nights of observations du ring the 2004 se a son. When the 2004 season data 
were combined with the observations of IWeldrake et al.l (120041 ) we were able to eliminate 
the proposed 5.25-d period, but the combined data set was insufficient to establish a unique 
ephemeris. We then examined the 286 V band images collected b y the OGLE team on 44 
nights during the 1993 observing season for the cluster field 104-A (IKaluzny et al.lll998l ). A 
single eclipse event was detected on the night of 1993 July 22 UT and this permitted the 
identification of an approximate orbital period of P ~ 29.540 d. Subsequent observations 
collected in the 2005, 2006, and 2007 seasons concentrated on nights with predicted eclipse 
events. During the 2007 season we observed the variable with the 2.5-m du Pont telescope 
using the 2048 x 2048 pixel TEK5 CCD camera at a scale of 0.259'7pixel. These data 
included one well covered primary eclipse observed mostly in the F-band on 2007 August 18 
UT. 

In addition to the Las Campanas photometry, we obtained V-band observations of an 
eclipse on 2005 October 22 UT using the 1.0-m telescope at the South African Astronomical 
Observatory. Observations were obtained with the STE4 1024 x 1024 pixel CCD camera at 
a scale of 0.31"/pixel. 

Prof ile photometry was extracted for all of the observations using the DAOPHOT/ALLSTAR 
package (IStetsonI 119871 ). We calibrated the instrumental photometry using observations col- 
lected with the du Pont telescope on the night of 200 7 August 19 U T. Observations were 
made of V69 together with 4 Landolt standard fields (jLandoltlll992l ). The standards were 
observed with a range of air-mass of 1.19 < X < 1.94. The conditions were photometric and 
the seeing ranged from 1.05 to 1.40 arcsec with a median value of 1.22 arcsec. The images of 
the cluster itself were obtained at an air-mass of 1.38 with seeing of 1.0 arcsec. Profile pho- 
tometry for the field of the variable was extracted from subframes covering 220 x 180 arcsec 
(the full field of TEK5 camera is 8.65 x 8.65 arcmin). This helped to minimize the effects 
of a variable point-spread function and as a result to obtain reliable aperture corrections. 
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Aperture corrections for the V69 obse rvations and the standard field observations were de- 
rived using the program DAOGROW (IStetsonlllQQOl ) . Magnitudes of 28 standard stars in the 
Landolt fields were taken from the Stetson catalog (Stetson 2000 1^. The following relations 
between instrumental (lower case letters) and standard magnitudes were obtained: 



v = V - 0.026(2) X (5 - 1/) + 0.118(3) x X + const, 
b = B- 0.069(2) x{B-V) + 0.217(3) x X + const, 



(1) 
(2) 



where X is the air-mass. In Fig. [T] we show the residuals between the standard and recovered 
magnitudes for the Landolt primary standards. The analysed V69 field includes 20 additional 
secondary standard stars from the Stetson catalog (IStetsonI 120001 ). The average residuals for 
these stars are AV = +0.004 ± 0.019 and A{B -V) = -0.001 ± 0.015 with our magnitudes 
fainter and our colors bluer on average. 

A total of 8 primary and 6 secondary eclipses were observed in the combined data sets. 
In Fig. [2] we show the BV light curves of V69 phased with the ephemeris : 



Mini = HJD 2453237.8421(2) + 29.53975(1) (3) 

The period was derived using the Lafier-Kinman algorithm, and the error in the period was 
estimated by visual inspection of the phased light curve as the period was varied away from 
the best fit value. These light curves contain a total of 1216 and 310 data points for V 
and B, respectively. The plots include photometry from all four data sets mentioned above. 
For the OGLE data we have plotted only points inside the primary eclipse. The colors and 
magnitudes of V69 at minima and at quadrature are listed in Table [H The quoted errors do 
not include possible systematic errors of the zero points of the photometric solution which 
we estimate to be 0.010 mag. 



SPECTROSCOPIC OBSERVATIONS 



We obtained spe ctroscopic observations of V69 with the MIKE echelle spectrograph 
(IBernstein et al.ll2003l ) on the Magellan Clay 6.5-m telescope. All observations were taken 
with a 0.7 arcsec slit at a resolution of R ~ 40,000. The observations generally consisted of 
two exposures fianking an exposure of a thorium-argon hollow-cathode lamp. Total exposure 
times per spectrum ranged from 1300 to 3665 seconds depending on observing conditions. 



^We have used the electronic version of the catalog as of 2007 Novem- 
ber 7. The catalog is maintained by Canadian Astronomy Data Centre at 



http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETSON/standards /I 
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The observa tions •were red uced with pipehne software written by Dan Kelson following the 
approach of iKelsonI (120031 ). Post-extraction processing of the spectra was done within the 
IRAF ECHELLE package! 



Velocities were measured with the TODCOR algorithm (jZucker fc Mazeh Ill994l ) using 
an implementation written by G. Torres. For velocity templates we used synthetic echelle 
resolution spectra from the grid of ICoelho et al.l (120061). We adopted the values of log g and 
Teff derived from the photometric solution (see next section) and assumed a metallicity of 
[Fe/H] = —0.71 with an a-element enhancement of 0.4 (see the discussion of the metallicity 
of 47 Tuc in the next section). The templates were adjusted once during the iterations 
for the Wilson-Devinney solutions, and final linear interpolations on the Coelho et al. grid 
were made at (log g, T^ff and [Fe/H]) = (4.14, 5945 K, —0.71) for the primary and (4.24, 
5955 K, —0.71) for the secondary. The measured velocities are insensitive to minor changes 
in these parameters. The interpolated synthetic spectra were smoothed with a gaussian 
to match the resolution of the observations. No rotational broadening was applied to the 
templates. The cross-correlations covered the wavelength intervals 4125 A < A < 4320 A, 
4350 A < A < 4600 A, and 4600 A < A < 4850 A. The final adopted velocities are the 
averages of these three measurements for each of the observations. 

The measured radial velocities were fit with a non-linear least squares solution using 
code written by G. Torres adopting the ephemeris given in equation (3). The observations 
are presented in Table [2] which lists the heliocentric Julian Date (HJD) at mid-exposure, 
the velocities of the primary and secondary components, and the orbital phases of the ob- 
servations. The adopted orbital elements are listed in Table |3] and the orbit is plotted in 
Fig. [31 



LIGHT CURVE ANALYSIS AND SYSTEM PARAMETERS 



We have used two different models for the aiialysis of the light c urves. The first 
is the Wilson-Devinney m odel (Wilson fc Devinnev 1971 : Wilson 1979 ^) as implemented 
in the PHOEBE package (jPrsa fc Zwittei 20051 ) . The secon d is the JKTEBOB program 



(jSouthworth et al. 



1981 



Etzel 



2004al lb[). which is based on the EBOP code JPoppeilEgSol : IPopper fc Etzel 



19811) . The most recent public version of JKTEBOB is described in detail in 



Southworth et al.l J2007h . In particular, it incorporates an option to adopt a non-linear limb 
darkening law. 



■^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 
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The two light curve fitting programs utilize different methodologies and approaches. 
The PHOEBE/WD program is more sophisticated in its treatement of the geometry and 
the stellar atmospheres. It utilizes Roche geometry to approximate the shapes of the stars, 
uses Kurucz model atmospheres, and treats the reflection effect in detail. It can be used 
for any component separation including contact systems. JKTEBOP/EBOP approximates 
stars by bi-axial ellipsoids. It is faster than the WD code but is appropriate only for binary 
components that are spherical or only slightly distorted. The real advantage of JKTEBOP is 
that the errors for the fit parameters are reliably determined. The PHOEBE code provides 
only formal errors of the fits and these are known to be underestimated (A. Prsa, personal 
communication). However an advantage of the PHOEBE code is that it can be used for the 
analysis of multi-color light curves. For V69 the solution for the 5-band is much poorer than 
for the \^-band because of the lower quality of the 5-band data. PHOEBE/WD was used 
to simultaneously fit the V and B curves. The resulting geometrical parameters are based 
mostly on the higher quality V curve while the simultaneous fit permits a good estimate of 
the luminosity ratio. 

In the analysis we assumed that the light curve of V69 is free from any "third light" . 
This is supported by the depths of the eclipses, with AV > 0.63 mag for both the primary 
and secondary eclipse. There is also no evidence for any third component in the spectra of 
the binary. Finally, we note that V69 can be located on numerous HST/ACS images (for 
example, PEP ID 9018, P.I. G. De Marchi). Examination of these images indicates that 
our ground-based photometry does not suffer from any blending problems from unresolved 
visual companions to the variable. 

Four recent high resolution studies of the metallicity of 47 Tuc suggest a value of 
[Fe/H] ~ —0.7. As part of a study to estab l ish a globular cluster metallicity scale based 
on measurements of Fe II lines, iKraft fc Ivand (120031) reanalysed ir i a un if orm way equivalent 
width s in a total of 8 giants me a sured by iBrown fc WallersteinI (Il992l ) , iNorris fc Da Costa 
fll995h . and lCarretta fc GrattonI Jl997h . The y obtained a rang e of -0 61 to -0.69 for [Fe/H]i 



and a range of -0.56 to -0.70 for [Fe/H]n. lAlves-Brito et al.l (120051 ) used observations of 



5 giants to obtain [Fe/H] y = -0.66 and [Fe/H jiT = -0.69 with an a-element enhancement 
of about 0.3 dex. iKoch fc McWilhamI J2008h found [Fe/H]i = -0.76 and [Fe/H]n = -0.82 
with an g-eler nent enhanceme r it of 0.4 dex based on a study of eight giants and one dwarf 
star. Finally, ICarretta et aP J2009h measured [Fe/H]i = -0.77 and [Fe/Hjn = -0.80 with 
an a-element enhancement of -1-0.39 from observations of 11 giants. For this analysis, we 
have adopted a value of [Fe/H] = -0.71 together with an a-element enhancement of -1-0.4. 
Sensitivities of age estimates to these assumptions are explored in Sections 6.1 and 6.2. 



The effective temperature of the primary is needed for the PHOEBE light curve solution 
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and for both components to estimate the luminosities of the stars from the derived radii. In 
the absence of detailed infrared eclipse photometry we estimated the effective temperature 
of the primary, Tp, from the dereddened {B — V) and {V — I) color indices. We used the 
observed color of the binary, (B — V) = 0.548 (se e Table [H) together with a reddening of 
E{B -V) = 0.04 jHarrislll996l ICratton et aUboOsI : IPercival et al]l2002h to derive {B - V)o 



= 0.508. For t he (V — I) color we note that V69 is star number 14065 in the catalog of 
Kaluznv et all Jl998h where {V - I) = 0.731±0 025. W e shift this observed color by -0.026 
to put the obser vation on the systern of [Stetson] (120001 ) (IPercival et al.ll2002l ) and adopt the 
reddening law of lSchlegel et al.l (Il998l ) to obtain E{V — I) = 0.06 and a final value of (V — I)o 
= 0.671. The similarity of the {B — V) colors in and out of eclipse (see Table [1]) suggests 
that the effective temperature of the primary derived from these colors is close to the true 
value. 

Visual inspection of the position of the binary in a color-magnitude diagram suggests 
that the components are beginning to evolve away from the main sequence (see Fig. H]), and 
so we adopt an initial gravity for the primary of log g = 4.1. 

These adopted colors and parameters lead to estimates for the effective temperature of 
the primary of 6070 K and 5900 K, 6050 K and 5905 K, and 596 K an d 5780 K for the (B-V) 
and (y — I ) colors and the cali b ration s of IWorthey &: Led (120061 ) . iRamirez &: Melendez 
( 120051 ) ■ and IVandenBerg &: ClemI (120031 ) . respectively. We adopt a linear average of these 
estimates, resulting in an effective temperature for the primary of 5945 ± 105 K, where the 
quoted error is the standard deviation of the six temperature estimates. 

The three color-temperature relations considered here do not explicitly a ccount for a- 



eleme nt enhancement. Based on the synthetic color-temperature employed by iDotter et al. 
( 120071 ) ■ we estimate that the difference in T^ff between [a/Fe]=0 and -1-0.4 is ~10 K for 
the [Fe/H], T^ff, and log g of the primary. We expect this is an upper limit because the 
color-temperature relations should include some implicit dependence on [a/Fe]. In any case, 
the [a/Fe]-color uncertainty is small compared to that which arises between the different 
color-temperature relations or between B — V and V — I. 



We note that iKoch fc McWilliamI (120081 ) adopted an effective temperature of 5750 K 
for the turnoff star that they studied based on excitation equilibrium for Fe I lines, and it is 
of interest to ask how accurate our photometric estimate of the effective temperature of the 
primary of V69 actually is. We estimate the uncertainty in the reddening to be 0.010 mag, 
and this value coupled with our estimated uncertainty in the photometric calibration sug- 
gests that the uncertainty in the derived unreddened {B — V) colors of the components 
of V69 to be ~0.014 mag. The {V — I) color is somewhat less accurate. For a range of 
metallicity of [Fe/H] from —0.61 to —0.77 (see above), an estimated range in the gravity of 
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0.2 dex and a uncertainty in [B — V)q of 0.014, the Worthey-Lee calibration suggests an 
uncertainty in the de r ived e ffecti ve temperature of approxini ately 80 K. The cahbrations of 



Ramirez fc Melendea (120051 ) . and I VandenBerg fc Clenj (120031 ) indicate similar uncertainties. 
Considering the scatter about T^ff in the various calibrations we conservatively conclude 
that we can estimate the effective temperature of the components of V69 to an accuracy 
of about 150 K. These three calibrations predict bolometric corrections spanning the range 
-0.11 to -0.09 and we adopt = BC^ = -0.10 ± 0.01. 



4.1. Models Using the PHOEBE Code 

The PHOEBE implementation of the Wilson-Devinney model fits the orbital inclination 
i, the gravitational potentials flp and fig, the effective temperature of the secondary Tg, the 
eccentricity e, the longitude of the periastron u and the relative luminosities Lg/Lp in B 
and V simultaneously. We iterated on the solution, calculating subsequent values for the 
{B — V)q colors for the primary and secondary from the observed {B — V)o for the system 
and the modeled values for {Ls/Lp)b and {Ls/Lp)v at each iteration. The mass ratio was set 
to the spectroscopic value of g = masss/massp = 0.984. The iterations converged in two 
cycles. The final values of the fit are listed in Table |H The values for the relative radii are 
derived directly from the non-dimensional potentials Qi and Q2, and the spectroscopic mass 
ratio q. The luminosity ratios and the relative radii are insensitive to changes in effective 
temperature: {Ls/Lp)B, {Ls/Lp)v, Ts, and Vp all change by less than 0.3 per cent for a 
±150 K change in effective temperature. 



4.2. Models Using the JKTEBOB Code 

The JKTEBOB code optimizes the sum of the relative radii r-p + r^, the ratio k = Vg/rp, 
the ratio of the surface brightness of the two stars J, the orbital eccentricity e, and the 
longitude of the periastron u. These two last parameters were included in the analysis 
by fitting e cos u and e sin u. The mass ratio was fixed at the spectroscopic value of 
q = 0.984, and the gravity brightening exponent was set to 0.32. Note that the values for 
these two parameters have a negligible effect on the model light curve as the stars are very 
well separated and both are practically spherical. We adopted a squ are-root law for the l imb 



darkening and adopted theoretical limb darkening coefficients from Ivan Hamme I (119931 ) for 



Teff = 5959 K (adopted from the PHOEBE code solution) and [Fe/H] = -0.71. 

Table contains the values of the fit parameters along with other relevant parameters 
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derived from the solution. The last two rows list the reduced chi-squared of the fit and the 
rms of the residuals. Separate solutions were derived for light curves in the V and B bands. 
The uncertaint ies of individual paramet ers were estimated using Monte Carlo simulations 
as described in ISouthworth et al.l (l2004bl ). Ten thousand simulations were run for each light 
curve. The solution based on the V light curve is much better constrained than the solution 
for the B band. The residuals of the fit obtained with the JKTEBOB are shown in Fig. [51 



4.3. Adopted Stellar Parameters 

A comparison of Tables H] and shows that the two light curve synthesis codes give very 
similar results. There is also good agreement between the fitted values of e and u and those 
derived from the spectroscopic observations. The remarkably small errors for the fit radii 
and relative luminosities might seem suspicious at first given that the light curves suggest 
that the eclipses of V69 are partial. However, the secondary eclipse is very close to being 
total. Examination of the synthetic light curves shows that only 0.4% of the surface of the 
secondary component remains visible at the center of the secondary eclipse. 

The magnitudes for the individual components of the binary are derived from the lumi- 
nosity ratios from the light curve solution and the observed out-of-eclipse photometry. We 
obtain Vp = 17.468±0.010, Bp = 18.019±0.010, K = 17.724±0.010 and 5, = 18.268±0.010 
where the errors are dominated by the contribution from zero point uncertainties of our pho- 
tometry, estimated at 0.010 mag. The colors of both components are practically identical 
with {B - V)p = 0.551 ± 0.014 and {B - V)., = 0.544 ± 0.014. The position of the binary on 
the 47 Tuc color - magnitude diagram is shown in Fig. |H The primary component has left 
the main-sequence of the cluster and is starting to ascend onto the subgiant branch. The 
secondary also has begun to evolve and is among the bluest stars at the cluster turnoff. 

The absolute parameters of V69 obtained from our spectroscopic and photometric anal- 
ysis are given in Table O Because the JKTEBOP code provides better error estimates than 
PHOEBE we have adopted the \^-band relative radii from Table O Note, however, that 
the relative radii derived with the PHOEBE models agree at a level of 0.2% with those ob- 
tained with the JKTEBOP models. The luminosities of the components follow directly from 
l oq(L/La,) = 2 loqiR /Rg,) +4:loQ(T,ff/ Ta,) where w e have adopted Rq = 6.9598 x 10^ km 



flBahcall et al.l 120051 ) . and Tq = 5777K (|Necke| ll986|) . 
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MEMBERSHIP OF V69 AND THE DISTANCE OF 47 TUC 



Before using the derived stellar parameters of the components of V69 to estimate the 
age of 47 Tuc, it is appropriate to consider the cluster membership of the star. There 
are four lines of evid ence available. First, although the velocity of 47 Tuc is low at Vrad = 
— 18. 7 ±0.5 km sec~^ ( IGebhardt et al.lll995l ). and thus does not provide a strong discriminant 



against projected field stars, the center-of-mass velocity of V69 (7 = —16.36 km sec~^) is 
consistent with this velocity. At the location of the variable - about 6 arcmin from the 
cluster center - the velocity dispersion of cluster stars is about 8 km sec~^. Second, Fig. H] 
shows that the individual components of V69 lie on the cluster main sequence. Third, 
the binary is a proper motion member of 47 Tuc. The absolute proper motion of 47 Tuc 
has been measured to b e fia = 5.64± 0.20 mas/yr and //^ = — 2.05± 0.20 mas/yr 
([Anderson fc Kingll2003l ). V69 is located in Field F of that study. This field has subsequently 
become a standard calibration field for HST, and V69 is present in numerous ACS/WFC 
frames taken between 2002 and 2006. A quick study of the ~30s F606W exposures shows that 
the proper motion of V69 relative to the bulk of the cluster stars is 0.2 mas/yr in RA and -0.1 
mas/yr in Dec, a typical internal motion for a cluster member at this radius (J. Anderson, 
personal communication) . Finally, based on the absolute magnitudes listed in Table [6] and 
the apparent magnitudes listed in Sec. 14.31 we derive distance moduli for the components of 
V69 of (m — M)y = 13.35±0.11 and {m — M)v = 13.34±0.12 for the primary and secondary, 
respectively. The average of these values is {m — M)v = 13.35 ±0.08. We compa r e this value 
to some recent distance determinations for 47 Tuc in Table [3 While iBono et al.l (120081 ) have 
shown that there appears to be some systematic errors in the measurement of the distance 
to 47 Tuc by different techniques, y 69 clearly lies at the distance of 47 Tuc. We note in 
particular that iKaluzny et al.l (l2007bl ) find a distance modulus of (m — M)v = 13.40 ± 0.08 
for OGLE-228, another eclipsing binary in 47 Tuc. We conclude that V69 is a member of 
the globular cluster 47 Tuc. 



THE AGE OF 47 TUC 



The potential of using observations of eclipsing binar ies for a robust determination of 
globular cluster ages was advocated by iPaczyhskil (119971 ). In particular, he argued that 
the mass-luminosity-age relation be used rather than the mass-radius- age relation since the 
latter can be affected by inaccuracies in the stellar models related to the treatment of sub- 
photospheric convection. In the following we compare the properties of the components of 
V69 with five different sets of stellar evolution models. We have attempted to minimize the 
effects of metallicity and a-element enhancement by comparing models with uniform values 
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of [Fe/H] = -0.71 and, where possible, [a/Fe] = +0.4. 



Victoria- Regina models (jVandenBerg et al.l 120061 ). and the Yonsei-Yale models (IKim et al. 



In our analysis we use the Dartmout h Stellar Evolution D atabase (IDotter et al.ll2007h . 
the Padova models ( Salasnich et al. 200oh. the Te ramo models ( Pietrinferni et al. 2006h. the 



Table [8] compares the parameters of each m odel set used in the deriv ation of the ages 
of V 69. Entries for the s o lar ab undance scale refer to lGrevesse fc Noelsl (119931 ) (Z/X ~ 0.024)) 
and iGrevesse fc Sauvall (119981 ) (Z/X ~ 0.023). We note that solar models based on recent 
measurements of the solar abundance scale remain in conflict with helioseismology data (see 
Serenelli et al.l (120091 ) for a discussion). We compare the derived values of mass, luminosity, 
and radius for the components of V69 with the isochrones. For the Dartmouth, Victoria- 
Regina, and Yonsei-Yale models we used interpolation software provided by the different 
groups to generate the plotted isochrones. For the Padova models we linearly interpolated 
on log (age) to generate isochrones spaced by 0.5 Gyr for ease of comparison with the other 
models. 

We plot the components of V69 on the mass-luminosity and mass-radius planes for the 
Dartmouth, Padova, Teramo, Victoria-Regina, and Yonsei-Yale models in Fig. El Fig. [TJ 
Fig. [HI Fig. [9l and Fig. [TOl respectively. In each Figure isochrones are shown with a step 
of 0.5 Gyr. The measured ages for the components of V69 are summarized in Table [91 To 
calculate the random errors in the age estimates we assumed that the measurement errors 
in mass, luminosity and radius are uncorrelated. This is a reasonable assumption for the 
errors in luminosity are dominated by the adopted error in T^ff and the errors in radius 
are dominated by the photometric solution rather than the error in A sin i as measured by 
the orbital solution. The errors in mass, luminosity, and radius lead to uncertainties in age 
through comparison with the sets of isochrones, and the final random errors listed in Table [9] 
are the quadrature sums of the mass-radius and mass-luminosity age uncertainties. The 
uncertainties in the masses (about 0.6%) and in the luminosities (about 10%) give similar 
contributions to the errors of the derived ages, while the random errors in the age estimated 
from mass-radius-age relations are completely dominated by the measurement errors in the 
masses of the components. We emphasize that the plotted uncertainty in luminosity arises 
directly from the ±150 K uncertainty in the effective temperatures of the components of V69. 
The radii are measured with a better relative accuracy than the luminosities and therefore 



''The models and supporting software are available at the following websites: Dart- 
mouth (http://stellar.dartmouth.edu/~models/), Padova (http://pleiadi.oapd.inaf.it/), 
Teramo ( http:/7l93.204.1.62/index.html[ ), Victoria-Regina (http:www.cadc- 
ccda.hia-iha.nrccnrc.gc.ca/cvo/community/VictoriaReginaModels/), and Yonsei-Yale 
( |http: / / www.astro.yale.edu/demarque / yyiso.htmlj ) . 
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the ages can be derived with a random error a factor of two smaller than the ages derived 
from the mass-luminosity relations. Table [9] also presents the weighted average of the ages 
measured for the primary and secondary components of the binary. 

Finally, since the components are very close to the cluster turnoff (see Fig. Hj), the 
age of the system can be estimated from turnoff mass-age relations. For each isochrone we 
estimated the turnoff mass as the mass at maximum log T^ff in that isochrone. Based on 
the position of the components in the CMD we conclude that the secondary is at the blue 
extreme of the CMD, and that the primary has begun to evolve onto the subgiant branch 
at a slightly cooler effective temperature. Our mass measurements are accurate enough to 
resolve the progression of mass along an isochrone in this part of the CMD. Fig. [TT] shows the 
resulting relations, and these derived ages are also summarized in Table O For completeness, 
ages are also given for the primary assuming it is at the cluster turnoff. The quoted errors 
in these ages are derived assuming only an error in the measurement of the masses, and are 
measured from the values of the ages at the ±lcr extrema of the measured masses. This age 
measurement is obviously fairly crude, but serves as a consistency check on the other ages 
estimates. 

The ages of V69 determined from these relations show some considerable spread from 
one model set to the next, ranging from 11.1 to 13.7 Gyr. While we only have age estimates 
from two stars, we can draw some general conclusions. The first is that the close agreement 
in the age estimates for the primary and secondary for all models suggests that the errors in 
luminosity are overestimated. Second, the ages measured by the three methods also agree to 
within the random errors within each model. This suggests that concerns about the model 
radii arising from uncertainties in the treatment of convection are overstated, and that the 
models predict internally consistent luminosities and radii at the measured masses. Future 
studies of clusters with more binaries covering a wider range of mass will be important in 
addressing this second issue. 

We next investigate whether or not there are systematic differences in the different model 
sets that might explain the wide range in measured ages. Many variables will affect the 
isochrones, among the most significant are adopted a-element enhancement and the relative 
abundance distribution over the a-elements, helium abundance, mixing length and convective 
overshoot, and the treatment of diffusion and gravitational settling. None of the models 
considered include the effects of rotation. Two of these variables ([Fe/H] and a-element 
enhancement) are measurable from spectroscopic data, while the others are parameters of 
the models. 

In the following subsections we explore the sensitivity of the measured ages to [Fe/H], 
[a/Fe], diffusion and gravitational settling, and helium abundance. 
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6.1. Sensitivity to Adopted [Fe/H] 

Figure [T2l shows the sensitivity of the measured ages to model [Fe/H] for the Dartmouth 
and Yonsei-Yale models sets for a range of ±0.05 dex around our adopted value of -0.71 for 
each of the three methods of measuring age. These values are interpolated in the models at 
a fixed [a/Fe] = +0.4. The plotted values are a weighted average of ages measured for the 
primary and secondary. The slopes (in terms of A Gyr per 0.1 dex) are given in Table [TOl 
In all three cases the slopes for the two model sets are very similar, and Table [10] also gives 
the average slopes. 



6.2. Sensitivity to a-Element Enhancement 

We have estimated the effect of differing a-element enhacement by measuring the ages 
of the components of V69 from Dartmouth and Yonsei-Yale isochrones calculated to have 
differing a-enhancement in the range 0.2 to 0.4 at a fixed [Fe/H] = -0.71. The results are 
shown in Fig. [131 where we plot the weighted average of the ages measured for the primary 
and secondary components against model [a/Fe]. As we found for the sensitivity to adopted 
model [Fe/H], the slopes (in terms of A Gyr per 0.1 dex) are similar for the two model sets. 
The individual values and the averages are listed in Table [TOl 



6.3. Helium Abundance and Heavy Element Diffusion 

At the adopted [Fe/H], the models span a modest but important range in helium abun- 
dance (y), from Y = 0.243 for the Victoria- Regina models to 0.2555 for the Dartmouth 
models. The age estimates from mass-radius, mass-luminosity, and turnoff mass-age rela- 
tions for all five model sets are plotted in the upper panel of Fig. [TH An obvious trend in 
age with Y can be seen in Fig. [TH as expected. In order to further explore the dependence 
of age on F, stellar evolu t ion tr acks were computed by one of us (Dotter) using the code 



described bv iDotter et all fl2007l ) for [Fe/H] = -0.71 and [a/Fe] = +0.4 at Y= 0.24, 0.255, 
0.27, 0.285, and 0.30 for the measured masses and ±1-0" uncertainties of the primary and 
secondary components of V69 (see Table [6]). Ages were read from the tracks at the measured 
luminosities and radii of the components and plotted as open symbols in Fig. [TH for all but 
the F=0.30 case. The plotted track-based ages are a weighted average of the ages of the 
primary and secondary components. These ages are listed in Table [TH The variation in 
measured ages with helium abundance for these tracks closely follows that seen for the five 
sets of model isochrones. 
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Some of the scatter seen in the upper panel of Fig. [M] comes from two sources. First, 
the Victoria- Regina models are calculated for [a/Fe]=+0.3, while the other four model sets 
use +0.4. To account for this we adjusted the Victoria-Regina models by +0.49 Gyr and 
+0.67 Gyr for ages determined from the mass-luminosity and mass-radius relations, respec- 
tively (see Fig. [13] and Table [TOll. Second, the Dartmouth and Yonsei-Yale models include the 
effects of helium and heavy element diffusion and gravitational settling while the Padova, 
Teramo, and Victoria-Regina r nodels do nqtl^ As a general rule, diffu sion will lower the 
ages determined from models (jSalaris et al.l l2000l : IChaboyer et al.ll200ll ). To estimate the 
importance of diffusion in the measured ages, evolutionary tracks were also calculated by 
Dotter for Y = 0.255 but without diffusion. Ages determined from these tracks are plotted 
in the upper panel of Fig. [TH It might be expected that the effect of diffusion will depend 
on adopted helium abundance, and in the absence of a detailed comparison of diffusion in 
multiple model sets we adopt as a rough estimate a correction of -0.7 Gyr as measured from 
the Dartmouth tracks with and without diffusion. This correction was applied to the ages 
determined from the Padova, Teramo, and Victoria-Reg ina models. We note that evidence 
for diffusion is seen in the globular cluster NGC 6397 (IKorn et al.l 120071 : iLind et al.l l2008l ) 
based on a systematic variation of measured metal abundanc es from the turnoff throug h 
the giant branch. No such evidence was detected in 47 Tuc by iKoch fc McWilliamI (l2008l ). 
although the 47 Tuc sample only included one turnoff star. 



The corrected ages are plotted in the bottom panel of Fig. [141 There are two conclusions 
to be drawn from this figure. First, the measured ages continue to be tightly correlated with 
Y, and in particular, remaining residuals from this trend are not strongly related to [a/Fe], 
distribution of the a-element enhancements, or the mixing length. This correlation means 
that absolute ages of globular clusters based on individual stars in these clusters will be 
crucially dependent on an accurate measurement of the helium abundance of the cluster. 
This might seem impossible given the accuracy of var ious methods of determining cluster 
helium abundances (see, for example, ISandquistI l2000l ) . However, we note that for model 
sets with Y < 0.25, the ages measured from mass-luminosity and mass-radius relations differ 
systematically, with mass-luminosity relations predicting higher ages. This effect is reversed 
for ages measured from Dartmouth tracks for Y = 0.285. If the models are sufficiently 
accurate, it should then be possible to constrain both age and Y by requiring that the 
mass- luminosity and mass-radius ages be consistent. 



^In fact, the situation is more complicated than this becaus e the Dartmouth mode ls employ a treatment 
of diffusion that is inhibited in the outer 0.01 Mq of the star (jChabover et al.ll200lh . The Teramo models 
employ diffusion in the calibration of their solar model but not in models used to generate their metal poor 
models and thus their predicted ages are indirectly affected by diffusion. 
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We investigate the issue of helium abundance further in Fig. [15] where we plot the 
measured radii and luminosities of the components of V69 together with Dartmouth tracks 
calculated for the masses of the primary and secondary for a range values of the helium 
abundance. The measured values of the radii and luminosities of both the primary and 
secondary agree to within the uncertainties for Yp = 0.267!!^q;o29 and Yg = 0.2711o;q3^ with 
an average value of F = 0.269too2i- This value can be compared with estimates of the 
helium abundance estimated from the pop ulation rat i o R = Nhb/ Rrgb measured from 
ground-based photometry (Y = .216l^-^|g: ISandquistI (120001 )) and HST photometry (Y = 
0.240 ± 0.015: ISalaris et al.l ( 20041 )). The agreement is reasonable given the uncertainties in 
the models used to calibrate the R method. There are at least two possible ways to improve 
upon our estimate. First, Figure 15 shows that our estimate is completely dominated by 
errors in the luminosities of the components of V69. Near-infrared eclipse photometry of V69 
can be used to improve upon distenace estimates, and hence luminostiy estimates. Second, 
discovery of additional similar systems to V69 in 47 Tuc will help improve improve upon 
distance, age, and hehum abundance measurements for this cluster. 



Casagrande et al.l (120071 ) derived bolometric luminosities and effective temperatures 



from multi-band photometry and the infrared flux method for a sample of 86 K-dwarfs. 
They compared these values to Padova stellar isochrones to estimate the He to metal en- 
richment ratio (AY/AZ) in the solar neighborhood. Although their primary interest was 
to measure AY/AZ, Casagrande et al. found that, in an absolute sense, the implied He 
content of the most metal-poor stars was as l ow as Y^O.l-in marked disagreement with 
canonical values for the primordial He content (ISalaris et al.l |2004J : ISpergel et al.l 120071 ) . In 
contrast, we find no indication that current stellar evolution models are unable to reproduce 
the mass-radius relationship of V69. 



6.4. The Absolute Age of 47 Tuc 

As mentioned above, a measurement of the absolute age of 47 Tuc will have random 
and systematic errors. Our best age estimates derived from the model isochrones have 
statistical errors of about 0.25 Gyr for mass-radius relations (see Table |9]). Outside of 
remaining systematic errors arising from different model parameters, the main systematic 
errors arise from the empirical estimates for metallicity and a-element enhancement. We 
estimate the errors in both of these measurements to be 0.10 dex, and following the results 
in Table [TOl we adopt a systematic error of 0.85 Gyr, the quadratic sum of the contributions 
from metallicity and a-element enhancement. This leads to an absolute age estimate of 
11.25±0.21(random)±0.85(systematic) Gyr from Dartmouth isochrones adopting Y = 0.255, 
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[Fe/H] = -0.71, and [a/H] = 0.4. Note that the age estimate from the Teramo isochrones 
is very similar when a correction for diffusion has been apphed (see Table |9] and Fig [T^ . 
as are the ages derived from mass-luminosity relations for both the Dartmouth and Teramo 
models. If we adopt a helium abundance of Y = 0.27 then the Dartmouth tracks indicate 
an age of 10.21±0.21±0.85 Gyr. 

This age estimate is cons istent with other recent estimates of the age of 47 Tuc. For 
example, iGratton et al.l (120031 ) find an age of 11.2 ±1.1 by fitting model isochrones with no 
diffusion to the clu ster color-magnitude di agram, and an age of 10.8 ± 1.1 Gyr for models 
includir ig diffusion. ISalaris fc Weisd (120021 ) derive an age of 10.7 ± 1.0, also from isochrone 
fitting. IZoccah et al.l J200lh obta in an age of 13 ± 2.5 Gyr by fitting the white-dwarf cooling 



Grundahl et al. I J2002h concluded that the age of 47 Tuc is "slightly below 12 Gyr'' 



sequence. 

based on analysis of Stromgren photometry and isochrone fitting. Our determination of the 
age of 47 Tuc based on observations of V69 has a small statistical error compared to that 
derived from isochrone fitting. 



7. DISCUSSION AND SUMMARY 

To the best of our knowledge, these measurements of the masses and radii of the com- 
ponents of V69 are the first such high accuracy (better than 1%) measurements to be made 
for Population II stars. The binary is a member of the globular cluster 47 Tuc and so the 
determination of its distance and age applies to the cluster as well. We obtained a distance 
modulus of (m — M)y = 13.35 ± 0.08. The main source of error in the distance estimate 
is the calibration dependent estimate of effective temperatures which we used to derive the 
bolometric luminosities for the component stars. 

A comparison of the measured masses, luminosities, and radii of the components to 
stellar evolution models suggests that the age of the system and hence the globular cluster 
47 Tuc can be measured to a statistical accuracy of about 0.25 Gyr. However, it is important 
to understand the assumptions that go into any one model. In particular, the derived 
ages are very sensitive to the adopted helium abundance. We derive an age for 47 Tuc of 
11.25±0.21±0.45 Gyr for a helium abundance Y = 0.255 using Dartmouth model isochrones. 
All models give similar ages when the effects of helium abundance are taken into account. 

Comparison of Dartmouth evolutionary tracks calculated for the measured masses of 
the primary and secondary indicate that the helium abundance can be measured to an 
accuracy of about 0.03 for each of the components. We estimate a helium abundance of 
Y = 0.269]'^q;o2i for 47 Tuc. The measured masses, radii, and luminosities of the components 
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of V69 are consistent with Dartmouth models assuming [Fe/H] = -0.71, [a/Fe] 
Y = 0.27. 



+ 0.4, and 



The radii of both stars are known with high accuracy, and it is therefore possible to 
obtain a more accurate and robust distance determination based on the surface brightness 
method (IBarnes fc Evanslll976l : lLacylll977l : [Thompson et al.ll200l[ ). The e mpirical calibration 



of surface brightness relations for dw arf and subgiant stars is improving (jPi Benedettdll998 



Kervella et al.ll2004j : lBuermannll2006l ). and it is reasonable to imagine that a distance accurate 
to a few per cent can be measured with accurate radii and {V — K) colors. We are in the 



proces s of collecting near IR eclipse profile photometry of both V69 and V228 (iKaluzny et al. 



2007bl ) in order to measure the distances to these two binary stars in this way. These data 
will improve the estimates of the bolometric luminosities of the components and lead to a 
more accurate measurement of the helium abundance and hence the absolute age of 47 Tuc. 

Finally, we note that contributions to the errors in the radii are dominated by the 
photometric solution. Given the large inclination, the errors in the masses are completely 
dominated by the orbital solution. An identical doubling of the existing set of radial ve- 
locity observations leads to a 33% improvement in the mass estimates, and a subsequent 
improvement in the age estimates. The system is relatively bright, and the prospects are 
good that a substantial improvement in the measured masses can be achieved with further 
radial velocity observations. 
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Appendix 



In their catalog of 47 Tuc variables IWeldrake et al.l (120041 ) included a potentially inter- 
esting detached binary named V39. The object has an orbital period of 4.6 d and is located 
about 0.30 mag to the red of the cluster main-sequence on the V/V — / plane. We have 



- 19 - 



obtained BV images of the V39 field using the du Pont telescope. These images show that 
V39 is a close visual pair of stars separated by 1.5 arcsec. The brighter component has 
V ^ 18.1 and B — V ^ 0.99 and is located on aV/B — V CMD among the Small Magellanic 
Cloud (SMC) asymptotic giant branch stars. The fainter component has V ~ 19.2 and 
B — V 0.08 and is a candidate SMC upper main sequence star. Given the orbital period 
of V39 we propose that it is the fainter component of the blend which is the eclipsing binary. 
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Table 1. BV Photometry of V69 at Minima and Quadrature 



Phase 


V 


B 


B-V 


Max 


16.836(1) 


17.384(2) 


0.548(2) 


Min I 


17.507(3) 


18.083(12) 


0.576(12) 


Min II 


17.461(3) 


18.011(14) 


0.550(14) 
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Table 2. Radial Velocity Observations of V69 



HJD Vp V, Phase 



(- 24()()()()()) 


km ri 


-1 


km H 


-1 






53183, 


.88089 


-57, 


.33 


24, 


.10 


0, 


.173 


53201, 


.93116 


21, 


.23 


-55, 


.92 


0, 


.784 


53206, 


.81454 


-5, 


,58 


-28, 


.40 


0, 


.950 


53206, 


.93425 


-7, 


.00 


-27, 


.15 


0, 


.954 


53210, 


.82680 


-41, 


.09 


8, 


.46 


0, 


.085 


53210, 


.93896 


-42, 


.27 


9, 


.98 


0, 


.089 


53271, 


.78367 


-53, 


.96 


20, 


,75 


0, 


,149 


53274, 


.77303 


-59, 


.04 


26, 


.63 


0, 


.250 


53280, 


.65021 


-23, 


.44 


-9, 


.68 


0, 


.449 


53281, 


.65012 


— 14, 


.50 


— 18, 


.34 


0, 


.483 


53282, 


.69872 


-7, 


,02 


-26, 


.32 


0, 


.519 


53521, 


.93280 


12, 


.54 


-45, 


.81 


0, 


.617 


53580, 


.91364 


11, 


.43 


-46, 


.41 


0, 


.614 


53581, 


.85305 


16, 


.04 


-50, 


.05 


0, 


.646 


53582, 


.85552 


19, 


.69 


-53, 


.69 


0, 


.680 


53584, 


.80006 


22, 


.23 


-56, 


.51 


0, 


.745 


53585, 


.79822 


21, 


.34 


-55, 


.77 


0, 


.779 


53631, 


.71181 


-50, 


.28 


16, 


.17 


0, 


.334 


53633, 


.74856 


-36, 


.21 


1, 


.48 


0, 


.402 


53634, 


.76492 


-27, 


.12 


-6. 


.54 


0, 


.437 


53889, 


.93336 


-37, 


.97 


6, 


.09 


0, 


.075 


53891, 


.93240 


-52, 


.13 


20, 


.77 


0, 


.143 


53892, 


.93208 


-56, 


.74 


24, 


.68 


0, 


.177 


53893, 


.93261 


-59, 


.42 


26, 


.95 


0, 


.210 


53935, 


.93349 


14, 


.75 


-48, 


.09 


0, 


.632 


53937, 


.92995 


21, 


.18 


-55, 


.15 


0, 


.700 


53938, 


.92894 


22, 


.39 


-56, 


.40 


0, 


.734 
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Tablc 3. Orbital Parameters for V69 



Parameter 


Value 


P (days) 


29.53975^ 


To (HJD-240 0000) 


2453237.8421^ 


7 {km s"^) 


-16.71± 0.05 


Kp {km s~^) 


41.03±0.13 


Kg {km s~^) 


41.86±0.09 


e 


0.0563±0.0011 


uj (deg) 


149.15±1.86 


Up {km s~^) 


0.50 


(7s {km s~^) 


0.33 


Derived quantities: 




A sin i {Rq) 


48.301±0.095 


Mp siTV^ i (Mq) 


0.8762±0.0048 


Ms sin^ i {Mq) 


0.8588±0.0060 



^Ephemeris adopted from photome- 
try. 



Table 4. Results of the Light Curve Analysis for V69 Obtained with the PHOEBE Code 



Parameter 



i (deg) 
Tp{K) 
Ts (K) 
e 

u (deg) 

{Ls/ Lp)y 

{Lg/ Lp)b 

rms 

{V) (mmag) 
(B) (mmag) 



89.771± 0.009 
5945 
5959±3 
0.0556±0.0002 
149.72±0.25 
0.02727±0.00005 
0.02408±0. 00008 
0.7870±0.0012 
0.7889±0.0015 
10.0 
14.5 
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Table 5. Results of the Light Curve Analysis for V69 Obtained with the JKTEBOP Code 



Parameter 



V 



B 



Adjusted quantities 

k = r-p/rs 
i (deg) 
J 
e 

u (deg) 

Other quantities 



Lg/ Lp 

arms (mmag) 



0.05166± 0.00009 
0.8836 ± 0.0072 
89.768± 0.012 
1.0146±0.0031 
0.0567±0.0008 
147.7±1.3 

0.02722±0.00009 
0.02405±0.00012 
0.792±0.011 
10.0 



0.05151±0. 00021 
0.8961±0.0199 
89.769±0.028 
1.0022±0.0091 
0.0586±0.0025 
145.2±3.4 

0.02717±0. 00026 
0.02435±0.00033 
0.802± 0.028 
14.5 
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Table 6. Absolute Parameters for V69 



Parameter Value 



A (Re) 




48.30±0.12 


M, (Mo) 




0.8762±0.0048 


M, (Mo) 




0.8588±0.0060 


Rp (Re) 




1.3148±0.0051 


Rs (Rq) 




1.1616±0.0062 


Tp(K) 




5945±150 


Ts (K) 




5959±150 






1.94 ±0.21 






1.53 ±0.17 


Mvp {mag) 




4.12 ±0.11 


Mvs {mag) 




4.38 ±0.12 


log 5rp(cm s" 




4.143 ±0.003 


log gs{cm s" 




4.242 ±0.003 
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Table 7. Comparison of Distance Determinations to 47 Tuc 



Author 


Method 


(m - M)v 


Gratton et al. f2003) 


Main Seq. fitting 


13.50±0.08 


Percival et al. f2002) 


Main Seq. fitting 


13.37±0.11 


Grundahl et al. f 20021 


Main Seq. fitting 


13.33 ±0.04 ±0.1 


Zoccali et al. f2001) 


white dwarfs 


13.27±0.14 


McLauehlin et al. f2006) 


kinematics 


13.15 ±0.13 


Salaris et al. f2007) 


IR luminosity of HB 


13.18 ±0.03 ±0.04 


Kaluznv et al. f2007b) 


Eclipsing binary 


13.40 ±0.07 


Bono et al. f2008) 


Tip RG branch 


13.32 ±0.09 


Bono et al. r2008) 


IR lum. of RR Lyrae stars 


13.47 ±0.11 


This paper 


Eclipsing binary 


13.35 ±0.08 
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Table 8. Parameters of Stellar Evolution Models 



Paranieter 


Dartnioutli 


Pado\-a. 


Teramo 


Victoria. 


Yonsei-Yale 


[Fe/H] 


-0.71 


-0.70 


-0.70 


-0.705 


-0.71 


[a/Fe] 


+0.40 


0.40 


+0.40 


+0.30 


+0.40 


X 


0.738010 


0.742000 


0.736000 


0.750250 


0.748031 


Y 


0.255000 


0.250000 


0.256000 


0.243000 


0.244646 


Z 


0.006490 


0.008000 


0.008000 


0.006750 


0.007323 


[0/Fe] 


0.40 


0.50 


0.50 


0.30 


0.40 


mixing length 


1.938 


1.680 


1.913 


1.890 


1.743 


(Z/X)so;ar 


0.023 


0.024 


0.024 


0.024 


0.024 



Table 9. Ages in Gyr for the Components of V69 Derived from Model Isochrones 



Method 


Component 


Dartmouth 


Padova 


Teramo 


Victoria- Regina 


Yonsei-Yale 


Mass - Luminosity 


primary 
secondary 

average 


11.11+0.54/-0.58 
11.12+0.80/-0.93 

11.11+0.45/-0.49 


12.39+0.55/-0.62 
12.46+0.80/-0.89 

12.41+0.45/-0.51 


11.93+0.55/-0.62 
12.04+0.76/-0.93 

11.97+0.45/-0.52 


13.62+0.57/-0.66 
13.90+0.68/-0.89 

13.74+0.44/-0.53 


13.02+0.54/-0.62 
13.16+0.72/-0.81 

13.07+0.43/-0.49 


Mass - Radius 


primary 

secondary 

average 


11.25±0.26 
11.25±0.37 

11.25±0.21 


12.26±0.29 
12.34±0.40 
12.29±0.24 


11.97±0.31 
12.16±0.42 

12.04±0.25 


13.32±0.33 
13.54±0.46 

13.39±0.27 


12.63±0.29 
12.70±0.39 

12.65±0.23 


Turnoff Mass 


primary 
secondary 


10.91±0.24 
11.75±0.31 


11.81±0.27 
12.79±0.37 


11.19±0.25 
12.15±0.35 


12.50±0.29 
13.58±0.40 


11.57±0.26 
12.52±0.34 
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Table 10. Effect of [Fe/H] and [a/Fe] on Measured Age 



Method 


Model 


[Fe/H] 


[a/Fe] 






AGyr / 0.1 dex 


AGyr / 0.1 dex 


Mass - Luminosity 


Dartmouth 


1.09 


0.67 




Yonsei-Yale 


1.03 


0.67 




average 


1.06 


0.67 


Mass - Radius 


Dartmouth 


0.73 


0.44 




Yonsci-Yale 


0.79 


0.54 




average 


0.76 


0.49 


Turnoff Mass 


Dartmouth 


0.68 


0.43 




Yonsei-Yale 


0.78 


0.37 




average 


0.73 


0.40 



Table 11. Ages in Gyr for the Components of V69 Derived from Dartmouth Tracks 



Method 


Component 


Y = 0.24 


Y = 0.255 


Y = 0.255^ 


Y = 0.27 


Y = 0.285 


Luminosity 


primary 


13.51+0.64/-0.66 


11.77+0.56/-0.65 


12.27+0.56/-0.65 


10.20+0.52/-0.63 


8.77+0.51/-0.62 




secondary 


13.71+0.76/-0.87 


11.85+0.72/-0.86 


12.40+0.76/-0.93 


10.16+0.71/-0.90 


8.60+0.74/-0.93 




average 


13.59+0.49/-0.53 


11.80+0.44/-0.52 


12.32+0.45/-0.53 


10.19+0.42/-0.52 


8.72+0.42/-0.52 


Radius 


primary 


12.97±0.32 


11.53±0.28 


12.32±0.30 


10.23±0.25 


9.10±0.22 




secondary 


12.98±0.45 


11.48±0.40 


12.49±0.43 


10.13±0.36 


8.94±0.32 






;i2.97±0.2(i 


ii.oiio/i:-! 


12..'-')<S±0.2.1 


10.21±0.21 


9.0r)±0.18 



^'Tracks calculated without the effects of diffusion. 
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Fig. 1. — Plot of the color and magnitude residuals for the standard stars observed on the 
night of 2007 August 19. 
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Fig. 2. — The phased BV hght curves of V69. 
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Fig. 3. — Radial velocity observations of V69. Filled symbols represent data for the primary 
and open symbols are for the secondary. 
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Fig. 4. — Position of V69 in the BV CMD for 47 Tuc. The square gives the position for the 
combined hght and the triangles show the positions of each of the individual components. 
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Fig. 5. — The residuals of the fits to the hght curve obtained with JKTEBOP for the 
primary (bottom) and secondary echpse (top). The residuals for B are offset by 0.15 mag 
for clarity. 
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0.85 0.86 0.87 0.88 

Mass (Mq) 

Fig. 6. — The masses, luminosities, and radii of the components of V69 are compared to 
isochrones based on the Dartmouth models. For each panel the isochrones are plotted in 
steps of 0.5 Gyr, with the lowest and highest age isochrones labeled. 



-41 - 




0.85 0.86 0.87 0.88 

Mass (Mq) 

Fig. 7. — The masses, luminosities, and radii of the components of V69 are compared to 
isochrones based on the Padova models. For each panel the isochrones are plotted in steps 
of 0.5 Gyr, with the lowest and highest age isochrones labeled. 
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Mass (Mq) 

Fig. 8. — The masses, luminosities and radii of the components of V69 are compared to 
isochrones based on the Tcramo models. For each panel the isochrones are plotted in steps 
of 0.5 Gyr, with the lowest and highest age isochrones labeled. 
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Mass (Mq) 

Fig. 9. — The masses, luminosities, and radii of the components of V69 are compared to 
isochrones based on the Victoria-Regina models. For each panel the isochrones are plotted 
in steps of 0.5 Gyr, with the lowest and highest age isochrones labeled. 



-44- 




0.85 0.86 0.87 0.88 
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Fig. 10. — The masses, luminosities, and radii of the components of V69 are compared to 
isochrones based on the Yonsei-Yale models. For each panel the isochrones are plotted in 
steps of 0.5 Gyr, with the lowest and highest age isochrones labeled. 
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Fig. 11. — Age as a function of turnoff mass for the Dartmouth (sohd hne), Teramo (dotted 
hne), Yonsei-Yale (long dashed hne), Padova (dot-dashed hne), and Victoria- Regina (short 
dashed hne) models. The masses for the primary and secondary components of V69 are 
plotted as solid vertical lines, with one-sigma errors represented by the shaded areas about 
the component masses. 



-46- 




Fig. 12. — Age as a function of [Fe/H] for Dartmouth (filled symbols) and Yonsei-Yale (open 
symbols) models. Ages are determined from mass- luminosity relations (bottom panel), mass- 
radius relations (middle panel), and turnoff mass - age relations (upper panel). 
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Fig. 13. — Age as a function of a-element enhancement for Dartmouth (filled symbols) and 
Yonsei-Yale (open symbols) models. Ages are determined from mass-luminosity relations 
(bottom panel) , mass-radius relations (middle panel) , and turnoff mass - age relations (upper 
panel) . 
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Fig. 14. — Age as a function of helium abundance as measured from Victoria-Regina, Yonsei- 
Yale, Padova, Dartmouth, and Teramo isochrones (sohd symbols plotted left to right, upper 
panel). Ages are determined from mass-luminosity relations (circles), mass-radius relations 
(squares), and turnoff mass - age relations (triangles). Open symbols represent ages measured 
from Dartmouth tracks, see text for details. Ages from two sets of Dartmouth tracks are 
plotted for Y = 0.255, models including hehum and heavy element diffusion (lower ages) 
and models without diffusion (higher ages). The plotted points have been offset in hehum 
abundance for clarity, the circles represent the correct helium abundance for any one model 
set. Bottom panel: Ages corrected for the effects of diffusion (Victoria-Regina, Padova, and 
Teramo models). The Victoria-Regina ages have been additionally corrected to bring the 
a-element enhancement to -1-0.4. 
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Fig. 15. — Dartmouth tracks calculated for the measured masses of the primary (left panel) 
and secondary (right panel). In each case the tracks are calculated for Y = 0.24 through 
Y = 0.30. The measured one sigma limits on mass are plotted as dotted lines for the Y = 0.27 
tracks. The measured values of radius and luminosity for the components of V69 are both 
consistent with Y = 0.27 with a one sigma range of 0.03. 



